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Abstract 

Two-dimensional almost-Riemannian structures are generalized Riemannian structures on 
surfaces for which a local orthonormal frame is given by a Lie bracket generating pair of vector 
fields that can become coUinear. We study the relation between the topological invariants 
of an almost-Riemannian structure on a compact oriented surface and the rank-two vector 
bundle over the surface which defines the structure. We analyse the generic case including the 
presence of tangency points, i.e. points where two generators of the distribution and their Lie 
bracket are linearly dependent. The main result of the paper provides a classification of oriented 
almost-Riemannian structures on compact oriented surfaces in terms of the Euler number of the 
vector bundle corresponding to the structure. Moreover, we present a Gauss-Bonnet formula 
for almost-Riemannian structures with tangency points. 

1 Introduction 

Let M be a two-dimensional smooth manifold. A Riemannian distance on M can be seen as the 
minimum-time function of an optimal control problem where admissible velocities are vectors of 
norm one. The control problem can be written locally as 

q = uX{q)+vY{q), + < 1 , (1) 

by fixing an orthonormal frame {X, Y). 

Almost-Riemannian structures generalize Riemannian ones by allowing X and Y to be collinear 
at some points. If the pair {X,Y) is Lie bracket generating, i.e., if 

span{X(<?), y (g) , [X, Y] (q) , [X, [X, Y]]{q),...}= T^M 

at every q E M, then (1) is completely controllable and the minimum-time function still defines 
a continuous distance d on M. Notice that a Riemannian distance can be globally defined on 
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M by a control problem as in (1) only if the Riemannian structure admits a global orthonormal 
frame, implying that M is parallelizable. More in general, X and Y are parallel on a set Z <Z M 
(called singular locus) , which is generically a one-dimensional embedded submanifold of M (possibly 
disconnected) . 

Metric structures that can be defined locally by a pair of vector fields {X, Y) through (1) 
are called almost- Riemannian structures. An almost-Riemannian structure S can be seen as an 
Euclidean bundle E of rank two over M (i.e. a vector bundle whose fibre is equipped with a 
smoothly- varying scalar product (•, and a morphism of vector bundles f : E ^ TM such that 
f{Eq) ^ TqM and the evaluation at q the Lie algebra generated by 

{f o a \ a section of E} 

is equal to TqM for every q G M. 

If E is orientable, we say that S is orientable. If E is isomorphic to the trivial bundle M x M^, 
we say that the almost-Riemannian structure is trivializable. The singular locus Z is the set of 
points of M at which f{Eq) is one-dimensional. An almost-Riemannian structure is Riemannian 
if and only ii Z = f/i, i.e. / is an isomorphism of vector bundles. 

The first example of genuinely almost-Riemannian structure is provided by the Grushin plane, 
which is the almost-Riemannian structure on M = with E = M? x M?, f{{x,y),{a,b)) = 
{{x,y), {a,bx)) and (•,•) the canonical Euclidean structure on M^. The model was originally in- 
troduced in the context of hypoelliptic operator theory [10, 11] (see also [4, 8]). Notice that the 
singular locus is indeed nonempty, being equal to the y-axis. Another example of (trivializable) 
almost-Riemannian structure appeared in problems of control of quantum mechanical systems (see 
[6, 7]). 

The notion of almost-Riemannian structure was introduced in [1]. In that paper, an almost- 
Riemannian structure is defined as a locally finitely generated Lie bracket generating C°°(M)- 
submodule A of Vec(M), the space of smooth vector fields on M, endowed with a bilinear, sym- 
metric map G : A X A ^ C°°{M) which is positive definite (in a suitable sense). This definition is 
equivalent to the one given above in terms of Euclidean bundles of rank two over M. A pair of vector 
fields {X, Y) in A is said to be an orthonormal frame for G on some open set Q if G{X, Y)(q) = 
and G{X, X){q) = G{Y,Y){q) = 1 for every g G 0. Equivalently, there exists a local orthonormal 
frame (a, p) for (•, •) such that X = foa,Y = fop. 

Almost-Riemannian structures present very interesting phenomena. For instance, even in the 
case where the Gaussian curvature is everywhere negative (where it is defined, i.e., on M \ Z), 
geodesies may have conjugate points. This happens for instance on the Grushin plane (see [1] 
and also [5] in the case of surfaces of revolution.) Moreover it is possible to define non-orientable 
almost-Riemannian structures on orientable manifolds and orientable almost-Riemannian structures 
on non-orientable manifolds (see [1]). 

This paper is a continuation of [1], where we provided a characterization of generic almost- 
Riemannian structures by means of local normal forms, and we proved a generalization of the 
Gauss-Bonnet formula. (For generalizations of Gauss-Bonnet formula in related contexts, see 
[2, 15].) Let us briefly recall these results. 

The flag of a submodule A of Vec(M) is the sequence of sub modules A = Ai C A2 C • • • C 
Am C • • • defined through the recursive formula 

Afc+i = Afc + [A,Afe]. 
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Denote by Am(g) the set {V{q) \ V £ Am}- Under generic assumptions, the singular locus Z has 
the following properties: (i) Z is an embedded one-dimensional submanifold of M; (ii) the points 
q € M at which A2((?) is one-dimensional are isolated; (iii) A3(g) = TgM for every q G M. We 
say that S satisfies (HO) if properties (i),(ii),(iii) hold true. If this is the case, a point q of M 
is called ordinary if A(q) = TgM, Grushin point if A(g) is one-dimensional and A2{q) = TgM, 
i.e. the distribution is transversal to Z, and tangency point if A2{q) is one-dimensional, i.e. the 
distribution is tangent to Z. If (fi, X, Y) is a local generator of A such that Q \ Z has exactly 
two components, then {X, Y) has different orientations on each of them. Normal forms for local 
generators of A around ordinary, Grushin and tangency points are recalled in Section 2.1. 

The main result of [1] is an extension of the Gauss-Bonnet theorem for orientable almost- 
Riemannian structures on orientable manifolds, under the hypothesis that there are not tangency 
points. More precisely, denote hy K : M\Z the Gaussian curvature and by u; a volume form 
for the Euclidean structure on E. Let dAg be the two-form on M \ Z given by the pushforward 
of uj along /. Fix an orientation S on M and let (respectively, M~) be the subset of M \ Z 
where dAs is a positive (respectively, negative) multiple of H. The main goal of [1] was to prove 
the existence and to assign a value to the limit 



where Mg = {q G M \ d{q,Z) > e} and d{-,-) is the distance globally defined by the almost- 
Riemannian structure on M. If M has no tangency point, then the limit in (2) can be shown to 
exist and to be equal to 27r(x(M"'') — x(M~)), where x denotes the Euler characteristic. When the 
almost- Riemannian structure is trivializable, we proved that xi^~^) = xi^~) whence the limit 
in (2) is equal to zero. Once applied to the special subclass of Riemannian structures, such a result 
simply states that the integral of the curvature of a parallelizable compact oriented surface (i.e., 
the torus) is equal to zero. In a sense, in the standard Riemannian construction the topology of 
the surface gives a constraint on the total curvature through the Gauss-Bonnet formula, whereas 
for an almost-Riemannian structure induced by a single pair of vector fields the total curvature is 
equal to zero and the topology of the manifold constrains the metric to be singular on a suitable 
set. 

The main objective of this paper is to complete the analysis in the more complicated case in 
which S has tangency points. 

The following result provides a classification of almost-Riemannian structures in terms of the 
Euler number of the vector bundle associated with it (for a definition of the Euler number see 
Section 2). 

Theorem 1 Let M be a compact oriented two-dimensional manifold endowed with an oriented 
almost-Riemannian structure S = (E, /, (•, •)) satisfying the generic hypothesis (HO). Then x(M+) — 
x{M~) + t{S) = e{E), where e{E) denotes the Euler number of E and t{S) is the number of rev- 
olutions of A on Z computed with respect to the orientation induced by M on Z. 

For a detailed definition of r(5) see Section 3. Notice that the Euler number e{E) measures how far 
the vector bundle E is from the trivial one. Indeed, E is trivial if and only if e{E) = 0. As a direct 
consequence of Theorem 1 we get that S is trivializable if and only if xiM^) — xiM~) -\- t{S) = 0. 




KdA 



(2) 
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For what concerns the notion of integrability of the curvature with respect to the Riemannian 
density on M \ Z, it turns out that the hypothesis made in [1] about the absence of tangency 
points is not just technical. Indeed, in Section 5.1 we present some numerical simulations strongly 
hinting that the limit in (2) diverges, in general, if tangency points are present. One possible 
explanation of this fact is the interaction between different orders in the asymptotic expansion 
of the almost-Riemannian distance. In this paper, to avoid this interference, we define a 3-scale 
integral of the curvature. Let T be the set of tangency points of S. Associate with every q £ T, 
a two parameters "rectangular" neighborhood B^^ {6i and 62 playing the role of lengths of the 
sides of the rectangle) built as follows. We consider a smooth curve (—1,1) 3 s w{s) passing 
through the tangency point w{0) = q and transversal to the distribution at q. We then consider, 
for each s £ (—1,1), the geodesic 7^ (parameterized by arclength) such that 7^(0) = w{s) and 
minimizing locally the distance from {w{s) \ s G (—1,1)} (such geodesic exists and is unique due 
to the transversality of the curve w; for details see [1]). For 81,62 sufficiently small, the rectangle 
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Figure 1: The rectangular box B'^^ 



B^^ is the subset of M containing the tangency point q and having as boundary 

752([-(5i,5i]) U7[„52_52](5i) U7_52([-(5i,(5i]) U7[„52,52](-(5i) 
(see Figure 1). Let Mir^Si,S2 — \ [JqeT 82- ^ 3-scale S-integrable if 

lim lim lim / KdAg (3) 

S^^O 52^0 6^0 J ^j^ ^^ ^^ 

exists and is finite. In this case we denote such a limit by KdAg. The following result, proved in 
Section 5.3, is a generalization of the classical Gauss-Bonnet formula for Riemannian structures to 
generic oriented two-dimensional almost-Riemannian structures. It concludes the analysis started 
in [1] including the presence of tangency points. 
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Theorem 2 Let M he a compact oriented two-dimensional manifold. For an oriented almost- 
Riemannian structure on M satisfying the generic hypothesis (HO), K is 3-scale S-integrahle and 



Notice that the construction of ,5^ is not intrinsic since it depends on the choice of a man- 
ifold transversal to Z at each tangency point and on its parameterization. The result, however, 
is. An interesting question is whether a canonical way of choosing these manifolds and their pa- 
rameterizations exists. This is related to the problem of finding intrinsic local normal forms for 
almost-Riemannian structures. 

As a consequence of Theorems 1 and 2 we get the following corollary. 

Corollary 1 Let M he a compact oriented two-dimensional manifold without boundary. For an 
oriented almost-Riemannian structure S on M satisfying the generic hypothesis (HO) we have 
^jyj KdAs = if and only if S is trivializable. In particular, if S has not tangency points then 
Jj^j KdAs = if and only if S is trivializahle. 

These results complete the analysis of the relation between the integral of the curvature and the 
topology of the manifold for two-dimensional almost-Riemannian structures. 

The structure of the paper is the following. In Section 2 we introduce the notion of n-dimensional 
rank- varying sub-Riemannian structure and recall some definitions and results given in [1]. In 
Section 3 we define the number of revolutions of a one-dimensional distribution along an oriented 
closed curve. In Section 4 we prove Theorem 1. First we construct a section of the sphere bundle 
on a tubular neighborhood of the singular locus having each tangency point as singularity. Then, 
we extend the section to the complement in the manifold of a finite set and show the sum of the 
indeces at all the singularities to be equal to xi^^^) ~ x{^^~) + 't{'S). In Section 5, the relation 
between tangency points and integrability of the curvature with respect to the area form associated 
with an almost-Riemannian structure is discussed. In particular, we provide in Section 5.1 some 
numerical simulations strongly hinting that in presence of tangency points the integral J^^ KdAg 
does not converge as e tends to zero. This leads us to introduce in Section 5.2 the notion of 
3-scale 5-integrability. Thanks to a Gauss-Bonnet formula for almost-Riemannian surfaces with 
boundary given in [9], we compute in Section 5.3 the total curvature of a generic two-dimensional 
almost-Riemannian structure with tangency points, proving Theorem 2. 

2 Basic definitions 

Let M be a n-dimensional manifold. Throughout the paper, unless specified, manifolds are smooth 
(i.e. C°°) and without boundary; vector fields and differential forms are smooth. Given a vector 
bundle E over M, the set of smooth sections of E', denoted by r(i?), is naturally endowed with the 
structure of C°°(M)-module. In the case E = TM we denote T{E) by Vec(M). 

Given an oriented vector bundle of rank n over a compact connected oriented n-manifold M, 
the Euler number of E^ denoted by e(£'), is the self-intersection number of M in E, where M is 
identified with the zero section. To compute e{E), consider a smooth section a : M ^ E transverse 
to the zero section. Then, by definition. 




e{E)= i{p,a) 



p\a(p)=0 
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where i{p,cr) = 1, respectively —1, if dpa : TpM — > T^^^-^E preserves, respectively reverses, the 
orientation. Notice that if we reverse the orientation on M or on E then e{E) changes sign. 
Hence, the Euler number of an orientable vector bundle E is defined up to a sign, depending 
on the orientations of both E and M. Since reversing the orientation on M also reverses the 
orientation of TM, the Euler number of TM is defined unambiguously and is equal to x{^)^ the 
Euler characteristic of M. 

Remark 1 Assume that a G ^{E) has only isolated zeros, i.e. the set {p \ cr{p) = 0} is finite. 
If E is endowed with a smooth scalar product (•, •), we define a : M \ {p \ cr{p) = 0} ^ SE by 

a{q) = -^=^M^=. Then if alp) = 0, i(p, a) = i(p, a) is equal to the degree of the map dB S"""^ 
V (o-{(j),CT(g)> 

that associate with each q G dB the value ^{q), where B is a neighborhood of p diffeomorphic to an 
open ball in that does not contain any other zero of a. 

Notice that if i{p,a) ^ 0, the limit Y\m.q^pa{q) does not exist and, in this case, we say that a 
has a singularity at p. 

Definition 1 A {n,k) -rank-varying distribution on a n-dimensional manifold M is a pair {E,f) 
where E is a vector bundle of rank k over M and f : E ^ TM is a morphism of vector bundles 
satisfying (i) / induces the identity on M i.e. f{Eq) C TqAl for every point q G M; (ii) the map 
o" I— > / o (J from T{E) to Vec(M) is injective. 

Remark 2 Definition 1 is equivalent to the definition of {n,k)-rank varying distribution given in 
[1], namely a submodule A o/Vec(M) locally generated by k (and not less than k) vector fields. 
Given a {n,k) -rank-varying distribution {E,f), we define A = {/ o cj | cr £ T{E)}. Thanks to (i) 
A is a submodule o/Vec(M). The condition (ii) implies that A cannot be generated by less than k 
vector fields. In the following, we use either definition, depending on our convenience. 

Let {E,f) be a (n, A;)-rank-varying distribution, A = {/ o a | a G ^{E)} be its associated 
submodule and denote by A(q) the linear subspace {V{q) \ V G A} = f{Eq) C TqM. Let Lie(A) 
be the smallest Lie subalgebra of Vec(M) containing A and Lieq(A) = {V{q) \ V G Lie(A)} for 
every q G M. We say that {E,f) satisfies the Lie bracket generating condition if Lieq(A) = TqM 
for every q £ M. 

A property (P) defined for (n, A;)-rank-varying distributions is said to be generic if for every 
vector bundle E of rank k over M, (P) holds for every / in an open and dense subset of the set 
of morphisms of vector bundles from E to TM inducing the identity on M, endowed with the C°°- 
Whitney topology. E.g., generically, a (n, A;)-rank- varying distribution is Lie bracket generating, 
provided that k > 1. 

We say that a (n, fc)-rank-varying distribution {E, f) is orientable if E is orientable as a vector 
bundle. Similarly, {E,f) is trivializable if E is isomorphic to the trivial bundle M x M'^. In this 
case, by definition, we have e{E) = 0. 

A rank- varying sub-Riemannian structure is defined by requiring that E is an Euclidean bundle. 

Definition 2 A {n,k) -rank-varying sub-Riemannian structure is a triple S = {E,f, {■,■)) where 
{E,f) is a Lie bracket generating {n,k) -rank-varying distribution on a manifold M and (•, •)<j is a 
scalar product on Eq smoothly depending on q. 

Remark 3 Definition 2 is equivalent to Definition 4 in [1], with G : A x A ^ C°°{M) defined by 
G{V,W) = {av,crw) where ay, aw are the unique sections of E satisfying /ouy = V,foaw = W. 
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Definition 2 generalizes several classical structures. First of all, a Riemannian structure {M,g) 
is a (n, n)-rank- varying sub-Riemannian structure with E = TM, f = Itm and (•,•) = g{-,-)- 
Classical sub-Riemannian structures (see [3, 14]) are (n, A;)-rank-varying sub-Riemannian struc- 
tures such that E' is a proper Euclidean subbundle of TM and / is the inclusion. Finally, we 
define n-dimensional almost- Riemannian structures (n-ARSs for short) as (n, n)-rank-varying sub- 
Riemannian structures. 

Let S = {E, f, (•, •)) be a (n, A;)-rank-varying sub-Riemannian structure. For every q G M and 
every v £ A{q) define 

Gg{v) = mi{{u,u)g \u£ Eg, f{u) = v}. 

If fJi, CTfc is an orthonormal frame for (•, •) on an open subset 0, of M, an orthonormal frame 
for G on Q, is given by / o ci ,...,/ o Cfc . One easily proves that orthonormal frames are systems of 
local generators of A. Notice that S is trivializable if and only if there exists a global orthonormal 
frame for (•,•). Hence S is trivializable if and only if there exists a global orthonormal frame for G. 

A curve 7 : [0, T] ^ M is said to be admissible for S if it is Lipschitz continuous and there 
exists a measurable essentially bounded function 

[0,T] 3t^u{t) G 

called control function, such that ^{t) = f{u{t)) for almost every t G [0,T]. Given an admissible 
curve 7 : [0, T] — > M, the length of ^ is 

rp 

%) = ^ ^G^W(7(0) dt. 

The function £(7) is invariant under reparameterization of the curve 7. Moreover, if an admissible 
curve 7 minimizes the energy functional J{'~i) = Jq G^^^) (7(i)) dt (with fixed T and fixed endpoints) 

then V = ^ G^(t) (7(i)) is constant and 7 is also a minimizer of ^(•). On the other hand a minimizer 
7 of such that v is constant is a minimizer of J(-) with T = i{^)/v. A geodesic for S is an 
admissible curve 7 : [0, T] — > M such that for every sufficiently small interval [ti, ^2] C [0, T], "y\[ti,t2] 
is a minimizer of J(-). A geodesic for which G^(f)(7(t)) is (constantly) equal to one is said to be 
parameterized by arclength. 

The distance induced by 5 on M is defined as 

d{qo, qi) = inf{^(7) I 7(0) = go, l{T) = gi, 7 admissible}. (4) 

The finiteness and the continuity of d{-, •) with respect to the topology of M are guaranteed by the 
Lie bracket generating assumption on the rank- varying sub-Riemannian structure (see [12]). The 
distance d{-, •) endows M with the structure of metric space compatible with the topology of M as 
smooth manifold. 

2.1 Normal forms for generic 2-ARSs 

Given a 2-ARS S, its singular locus is the set of points q where A(g) = f{Eq) has not maximal 
rank, that is, 

Z = {qeM\ dim(A(g)) = 1}. 
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Notice that dim(A(q')) > for every point q S M, because A is a bracket generating distribution. 

We say that S satisfies condition (HO) if the fohowing properties hold: (i) Z is an embedded 
one-dimensional submanifold of M; (ii) the points g G M at which A2(g) is one-dimensional are 
isolated; (ill) A3(g) = T^M for every q £ M, where Ai = A and A^+i = A^ + [A, A^]. 

Proposition 1 ([1]) Property (HO) is generic for 2- ARSs. 

ARSs satisfying hypothesis (HO) admit the following local normal forms. 

Theorem 3 ([1]) Given a 2-ARS S satisfiyng (HO), for every point q € M there exist a neigh- 
borhood U of q and an orthonormal frame {X, Y) for G on U , such that up to a smooth change of 
coordinates defined on U, q = (0, 0) and {X, Y) has one of the forms 



(Fl) 


X{x,y) 


= (1,0), 


Y{x, 


.y) 


= (0,e'^(^'^)), 


(F2) 


X{x,y) 


= (1,0), 


Yix; 


~y) 


= (0,xe<^(^'^)), 


(F3) 


X{x,y) 


= (1,0), 


Y{x, 




= (0, {y - xV(a:))e«(^'J') 



where (p, ^ and ip are smooth real-valued functions such that 0(0, y) = and V'(O) > 0. 

Let 5 be a 2-ARS satisfying (HO). A point g G M is said to be an ordinary point if A(g) = TqM, 
hence, if S is locally described by (Fl). We call q a Grushin point if A{q) is one-dimensional and 
A.2{q) = TgM, i.e. if the local description (F2) applies. Finally, if A(g) = A2{q) has dimension one 
and A3(g) = TgM then we say that g is a tangency point and S can be described near q by the 
normal form (F3). We define 

T = {q £ Z \ q tangency point of S}. 
2.2 A Gauss— Bonnet formula for 2- ARSs 

Let M be an orientable two-dimensional manifold and S = {E, f, {■, •)) an orientable 2-ARS on M. 
Notice that (•, •) defines a Riemannian structure on M \Z. Denote by K the Gaussian curvature 
of such a structure and by u; a volume form for the Euclidean structure on E. Let dAg be the 
two-form on M \ Z given by the pushforward of uj along /. Once an orientation on M is fixed, 
M \ Z IS split into two open sets M"*" and M~ such that the orientation on M coincides with the 
one defined by dAs on M'^ and with its opposite on M~. 

For every e > let = {q G M \ d{q, Z) > e}, where d{-, •) is the almost-Riemannian distance 
(see equation (4)). We say that K is S-integrable if 

lim / K dAs 

exists and is finite. In this case we denote such limit by Jj^j KdAg. 

When S has no tangency points K happens to be 5-integrable and Jj^^ KdAg is determined by 
the topology of M'^ and Af~. This result, recalled in Theorem 4, can be seen as a generalization 
of Gauss-Bonnet formula to ARSs. 

Theorem 4 ([1]) Let M he a compact oriented two-dimensional manifold, endowed with an ori- 
ented 2-ARS S for which condition (HO) holds true. Assume that S has no tangency points. Then 
K is S-integrable and Jj^.j KdAs = 27r(x(M+) - x{M~)). 
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3 Number of revolutions of A 



Let M be a compact oriented two-dimensional manifold and W an oriented closed simple curve in 
M. Since M is oriented, TM\w is isomorphic to the trivial bundle x and its projectivization 
is isomorphic to W x S^. Every subbundle T of TM\w of rank one can be seen as a section of the 
projectivization of TM\]y i.e. a smooth map T : W ^ W x such that vri o T = Idw, where 
TTi : W X ^ W denotes the projection on the first component. We define r(T, W), the number of 
revolutions of T along W, to be the degree of the map 7r2 o T : 1^ — > S^, where 112 W x ^ 
is the projection on the second component. Notice that t(T,W) changes sign if we reverse the 
orientation of W. 

Remark 4 Let M,W,T as above and let us show how to compute t{T,W). Let V G T(TW) be 
a never-vanishing vector field along W. Then span(y) is a subbundle ofTM\w of rank one and 
there exists an isomorphism t : TM\\y — > x such that t o V{q) = (g, (1, 0)). The trivialization 
t induces an isomorphism, still denoted by t, between the projectivization ofTM\\Y and W x 
such that 1:2 ° t o V : W ^ is constant. To simplify notations, we omit t in the following. Let 
■K2 o V{q) = 6q. Assume that Oq is a regular value o/7r2 o T (otherwise there exists a smooth section 
T homotopic to T having 6q as a regular value). By definition, 

t{T,W)= Y1 signK(vr2oT)), 

i}|7r2oT((j)=6»o 

where dq denotes the differential at q of a smooth map. Notice that a point q satisfies 7r2oT(g) = 6q 
if and only if T(q) is tangent to W at q. 

Let M be a compact oriented two-dimensional manifold and S be an oriented 2-ARS satisfying 
(HO). Then the singular locus Z is the boundary of M"'". Fix on Z the orientation induced by 
M"*" and let C{Z) denote the set of connected component of Z. Since A is one-dimensional along 
Z, we can define r(5) = Y^weCiZ) '^(^' 

Remark 5 Let S be a 2-ARS satisfying hypothesis (HO). Let V,6o,7r2 be as in Remark 4 with 
T = A and W £ C{Z). Since ^2,{q) = TqM for every point q G M , 9q is a regular value 0/712 o A. 
Moreover, the set of points q £ W such that 1:2 ° A((/) = 9q is the set of tangency points of S 
belonging to W . Hence, 

t{A,W)= Y1 sign(d,(7r2 o A)). 
qewnr 

We define Tq = sign(dg7r2 o A) (see Figure 2). Clearly, t{S) = Ylq^r'^i- 

4 Proof of Theorem 1 

The idea of the proof is to find a section a of BE with isolated singularities pi , . . . , pm such that 
X^jLi ^(Pii — x{^^) ~ x{M^) + t{S). In the sequel, we consider Z to be oriented with the 
orientation induced by M^. 

We start by defining cj on a neighborhood of Z. Let W he a, connected component of Z. 
Since M is oriented, there exists an open tubular neighborhood W of W and a diffeomorphism 
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A 


A 


/ / \ \ ^+ 

W / / . \ \ w 

/ / , \ \ M- 

r, = -1 


/ 

^ ' \ / / ^"^ 
\ W / M- 

Q 

r, = 1 



Figure 2: Tangency points with opposite contributions 



^' : 5^ X (—1,1) — > W that preserves the orientation and ^'l^ixjo} orientation-preserving 
difFeomorphism between and W. Remark that / : E\t^+ TW^ is an orientation-preserving 
isomorphism of vector bundles, while / : E\-^- TW" is an orientation-reversing isomorphism 
of vector bundles, where = WnM^. For every s 0, lift the tangent vector to 9 ^ ^(^j s) to 
E using /~^, rotate it by the angle tt/2 and normalize it: a is defined as this unit vector (belonging 
to if s > 0, its opposite if s < 0. In other words, a : W \ — > SE is given by 

(j(g) = sign(g) (6»,s) = ^' (g), (5) 



^{f-^{^{0,s))J~\%{e,s))) 



where -R,r/2 denotes the rotation (with respect to the Euclidean structure) in E by angle 7r/2 in 
the counterclockwise sense. The following lemma shows that a can be extended to a continuous 
section from W \ T to SE. 

Lemma 1 a can he continuously extended to every point q \ T . 

Proof. Let q ^ W \ T , U he a, neighborhood of g in M and {x,y) be a system of coordinates 
on U centered at q such the almost-Riemannian structure has the form (F2) (see Theorem 3). 
Assume, moreover, that [/ is a trivializing neighborhood of both E and TM and the pair of vector 
fields {X, Y) is the image under / of a positively-oriented local orthonormal frame of E. Then 
W r\U = {{x^y) I X = 0}. Since ^(0,0) is non-zero and tangent to W, ^{9,0) is tangent to the 
y-axis. Hence, thanks to the Preparation Theorem [13], there exist /12 : M — > M, /ii,/i3 : — > M 
smooth functions such that /i2(y) 7^ for every y G M and for ^{9, s) £ U 

— (61,5) = {xhi{x,y),h2{y) + xhs{x,y)), 

where (x, y) are the coordinates of the point ^'(^, s). Let us compute a at a point p £ (WCiU) \ W. 
Since 

\ h f \vr \ , ^■iiv) +xh^{x,y) 
— {9, s) = xhiix, y)X{x, y) + Y{x, y), 

then 
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where {C, p) is the unique local orthonormal basis of E\u such that f o C, = X and f o p = Y . 
Notice that U D M"*" = {{x,y) | x > 0} and U D M~ = {{x,y) \ x < 0}. Using formula (5), for 
(x, y) = ^{9,s) € U\W one easily gets 

sign(x) / h2{y) + xh^{x,y) 
a(x,y] = —, r it — ^ (. + xrii(x,y)p 

where 



Since 



j.^ sign(x)(/i2(j/) + xh3{x,y)) ^ h2{y) ^.^ sign(x)x/ii(x, y) 



x^o l{x,y)xe'l'('^'y) I ^2(2/) I ^'^0 Kx,y) 

a can be continuously extended to the set {x = 0} = W HU. ■ 

The next step of the proof is to show that for every q £ W ClT , i{q, a) = Tq. 

Lemma 2 Let a : W \ T — > SE he the continuous section obtained in Lemma 1. Then, for every 
q £ W nT the index of a at q is equal to Tq and, consequently, 

i{q,a) = TiA,W). (6) 

q&rnw 

Proof. Let q £ W OT, U he a neighborhood of g in M and (x, y) be a system of coordinates 
on U centered at q such the almost-Riemannian structure has the form (F3) (see Theorem 3) i.e. 
a local orthonormal frame {X, Y) is given by 

X = {1,0), y = (0,(y-xV(x))e«(^''^)). 

Define a = 1, respectively —1, if (X,Y) is the image under / of a positively-oriented, respectively 
negatively-oriented, local orthonormal frame of E. One can check that Tq = —a. Let us make the 
following change of coordinates 

X = X, y = a{y — x^^(x)). 

In these new coordinates, X and Y become 

X = (1, -a(2x'0(5) + Y = (0, ye^i^'^v+^'^Pi^))) 

and W r\U IS the x-axis. In the following, to simplify notations, we omit the tildes and denote the 
function ^(x, ay + x'^ip{x))) by ^(x, y). Since §f 0) is tangent to W, by the Preparation Theorem 
[13] there exist hi : M — > R, /12, /13 : M'^ — > M smooth functions such that /ii(x) 7^ for every x £ M 
and for ^(9,3) £ U 

dip 

— (61, s) = (/ii(x) + yh2{x, y),yh3{x, y)) 
where (x,y) are the coordinates of the point ^'(6',s). This implies that 

dib 

-^(.0,s) = {hi{x)+yh2{x,y))X{x,y) + 

yh3{x,y) + a{hi{x)+yh2{x,y)){2xilj{x)+x'^\l)'{x)) 
+ Ti — \ Y[x,y). 
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Let (C, p) be the local orthonormal frame of E such that X = foC, and Y = fop. From equation (5), 
it follows that 

/ N sig-n.{y)yh:i{x,y) + a{hi{x)+yh2{x,y)){2xil}{x) +x^\Ij'{x)) 

a[x,y) = —a— 7 -, — r Q + 

^ ' l{x,y) ye^i^^'V) 

+ "Tr-^(^i(^) + yh2{x, y))p, 
l{x,y) 

where 



Kx, y) = \ {hi{x) + yh2{x, y)Y + 



Notice that for 2; = 0, y 7^ we have 



yhz{x,y) + a{hi{x) + yh2{x,y)){2xii{x) + x'^ip'{x)) 



ye 



^(0'^) = ..rlr^^wn ^2 2.( f^(-""^^°''^^^(°'^)^ + (hm+yh2{0,y))p), 

v/(/ii(0) + y/i2(0, y))2 + /l3(0, y)2e-2?(0,y) 

whence the limit of a as (x, y) tends to (0, 0) does not exist. Let us compute the index of a at 
q = (0,0). Using Taylor expansions of the components of a in the basis {C, p) we find 

Take a circle t 1-^ (rcosi, rsini) of radius r centered at (0,0) and assume r so small that (0,0) is 
the unique singularity of a on the closed disk of radius r. By definition, i((0, 0), a) is half the degree 
of the map from the circle to M/vrZ that associates to each point the angle between span(cr) and 
the C- Using (7), this angle is 

/ N f y/ii(0)e«(0'0) 

o(x, y) = — arctan — — 7 — — - — -— + Ci-v/a; + V ) • 

^ ^ \yhs{0,0) + 2ax^{0)hi{0) ^'j 

Computing a along the curve x{t) = rcost, y{t) = rsint, we find 

. ^ I sin(t)/ii(0)e«(0'0) \ 

air cos t, r smi = — arctan -— — ; ^ , , , , , , — -— + Uir] . 

^ ' ' \^sin(t)/i3(0,0) + 2acos(t)V'(0)/ii(0) ^'j 

Hence, by letting r go to zero, we are left to compute the degree of the map a : [0, 27r) [0, vr) 
where 

/ sin(t)/ii(0)e«(o.o) \ 

alt) = — arctan -— — ^ , , , , , , — —— . 

^ ' \^sin(t)/i3(0,0) +2acos(t)^(0)/ii(0) j 

Since zero is a regular value of a, the degree of a is 

sign(a'(t)) = sign(a'(0)) + sign(a'(7r)) = —2a, 

iG[0,27r)|a(t)=0 

where the last equality follows from a'(0) = a'(vr) = — ae^^''''')/(2'0(O)). Hence, i(g, o") = —a. Since 
Tq = —a, the lemma is proved (see Remark 5). ■ 
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Let C{Z) denote the set of connected component of Z. Let Z = \J\y^c{z) '^^ ^^'^ consider 
an orientation-preserving diffeomorphism ^ : Z x (—1,1) \lweC{Z)^ such that ^|jx{o} 
an orientation-preserving diffeomorphism onto Z. Applying Lemma 1 to every W £ C{Z) and 
reducing, if necessary, the cyhnders W, we can assume that the set of singularities of o" on U = 
L[vKeC(2) ^ is ^ ■ Then a : U \ T — > SE is continuous. Moreover, by equation (6), 

^i{q,(j) = t{S). 

Extend o" to M \ U. By a transversality argument, we can assume that the extended section has 
only isolated singularities {pi, . . . ,pk} £ M \ Z. Since 

k k 
j = l q£T j = l 

we are left to prove that 

fc 

Y,i{pj,a) = x{M-^)-x{M~). (8) 
i=i 

To this aim, consider the vector field F = f oa. F satisfies G{F, F) = 1, where G(-, •) is defined 
as in Remark 2 and the set of singularities of -F|m\2 is exactly {pi, . . . Let us compute the 

index of F at a singularity p E {pi, . . . ,Pk}- Since / : TM^ preserves the orientation and 

/ : E\]^f~ — > TM~ reverses the orientation, it follows that i{p,F) = ±i{p,a), ifp G M^. Therefore, 

k 

j=l j\PjeM+ j\vj&M- 

The theorem is proved if we show that 

i{pj,F) = x{M+), i{Pj,F)=x{M-). (10) 

To deduce equation (10), define = M+ \ ^{Z x (0,1/2)). Notice that, by construction, 
'^\^{Zx{i/2}) is non-singular, hence the same is true for -^|,i,(Jx{i/2})- Moreover, the almost- 
Riemannian angle between Tq{dN'^) and span(F(g)) is constantly equal to 7r/2. Hence F\Qf^+ 
points towards and applying the Hopf 's Index Formula to every connected component of A^"*" 
we conclude that 

j\pj£M+ j\Pj&N+ 

Similarly, we find 

Y i{Pj,F)=x{M-). 

j\pj&M- 
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5 <S-integrability in presence of tangency points 



5.1 Numerical simulations 



In this section we provide some numerical simulations hinting that, when T 7^ 0, 



/ KdA, 



■s 



J Me 



does not converge, in general, as e tends to zero. 

Prom the proof of Theorem 4 we know that far from tangency points the integral of the geodesic 
curvature along dM^ and dM~ offset each other for e going to zero. Hence, to understand whether 
the presence of a tangency point may lead to non-5-integrability of K it is sufficient to compute the 
geodesic curvature of dM^ and dM~ in a neighborhood of such a point. More precisely consider 
the almost-Riemannian structure {E^ /, (•, •)) on M = for which = x M'^, /((x, y), (o, h)) = 
((x, y), (a, h{y — x^))) and (•, •) is the canonical scalar product. For this system one has 



liminfq^(Q 0) -fir(q') = —00. This situation is different from the Grushin case where K{q) diverges 
to —00 as g approaches Z. 

For every e > 0, the sets dM^ and dM~ are smooth manifolds except at their intersections 
with the vertical axis x = 0, which is the cut locus for the problem of minimizing the distance 
from Z = {(x,x^) | x G M}. Fix < a < 1 and consider the two geodesies starting from the 
point (a, a^) and minimizing (locally) the distance from Z. Let P'^ and P~ be the two points 
along these geodesies at distance e from Z. Denote by 7"*" and 7" the portions of dM^ and 
dM~ connecting the vertical axis to the points and P~, oriented as in Figure 4. It is easy 
to approximate numerically 7"*" and 7" by broken lines, but the evaluation of the integral of their 
geodesic curvatures is very unstable since its computation involves the second derivative of the 
curve parameterized by arclength. To avoid this problem, we rather apply the Riemannian Gauss- 
Bonnet formula on the regions Q.'^ and introduced in Figure 4. This works better since the 
integral of the Gaussian curvature on O"*" and 0" is numerically stable, and the integral of the 
geodesic curvature on horizontal and vertical segments can be computed analytically (in particular 
it is always zero on horizontal segments). Figure 5 shows the value of 



for a = 0.1 and e varying in the interval [0.01,0.04]. The graph seems to converge as e tends to 
zero to a nonzero constant, strongly hinting at the divergence of Jj^^ KdAg. 




+00 and 



£ 




14 



Figure 3: Graph of K for A = span((l, 0), (0, y — x'^)) 
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Figure 5: Divergence of the 5- integral of K 



5.2 More general notion of 5-integrability 

The simulations of the previous section lead us to introduce the following alternative notion of 
integrability. 

Definition 3 (3-scale 5-integrability) Let q £ T and If^ be a neighborhood of q such that an 
orthonormal frame for G on is given by the normal form (F3). For 5i, (52 > sufficiently small 
the rectangle [— x [—82,62] is a subset of W denoted by For every e > 0, define 



exists, is finite and does not depend on the choice of the normal form. In this case we denote such 
limit by ^j^j KdAg . 

Remark 6 Notice that if T = 0, then the concepts of S -integrability and 3-scale S -integrability 
coincide. 

The order in which the limits are taken in (12) is important. Indeed, if the order is permuted, 
then the result given in Theorem 2 does not hold anymore. 

Recall that the normal form (F3) is not totally intrinsic, since the functions and cj) depend 
on the choice of a parametrized smooth curve passing through the tangency point and transversal 
to the distribution at the point. 

5.3 Proof of Theorem 2 

Let us recall the following Gauss-Bonnet-like formula for domains whose boundary is in a 
neighborhood of Z. 




(11) 



We say that K is 3-scale S-integrable if 




(12) 
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Theorem 5 ([9], Theorem 5.2) Let U be an open bounded connected subset of M such that i) U 
contains only ordinary and Grushin points, ii) dU is piecewise , iii) dU is in a neighborhood 
of Z, iv) dU is the union of the supports of a finite set of curves 7^, . . . ,7"* that are admissible 
for A and of finite length. 

Define = n U. Then the following limits exist and are finite 



KdAs := lim 

£ 



1 / KdAs, (13) 
Ju+uur 



idl 



au \JdUndUt JdUndU, 

where we interpret each integral ^qij^qu± kgda as the sum of the integrals along the portions of 
dU n dU^, plus the sum of the angles at the points of dU n dU^ where dU is not C^. Moreover, 
we have 

[ KdAs + [ kgdos = 27r(x(C/+) - x{U~))- (15) 

JU JdU 

Fix 5i and 82 in such a way that the rectangles Bj^ are pairwise disjoint and Z n dB^^ C 

[— (5i,5i] X {52}, for every g G T. By construction, dBl^ is admissible and has finite length for 

every (7 G T. Hence we can take Mg \ UggT 52 ^ ^ Theorem 5. As a consequence, we have 

hm I KdAs + ^ / kgdas = 2n{x{M+ \ [J Bl ,J - x(M- \ [j Bl ,J) 

= 27r{x{M+)-x{M-)) = 27T{e{E)-T{S)), 
where the last equality follows from Theorem 1. We are left to prove that, for a fixed q € T, 

lim lim / kf,da<, = —2ttt„ 

-0 J 91 



kgdcFs := lim | / kgda — j kgda j , (14) 



"1 >"2 



(see Remark 5). In order to prove it, let us work with the normal form (-F3) and assume that 
M"*" n is the set {y — x^ip{x) < 0} H f/"^, the proof for the opposite situation being analogous. On 
one hand, one can check that Tq = 1. On the other hand, the geodesic curvature along [— (5i, 61] x {82} 
and along [— x {—82} is zero, the two segments being the support of geodesies. Hence 

4 

/ kgdas = / kgdas + / kgdas + V] 

where the last term is the sum of the values of the angles of the box and is equal to —2tt. Indeed, 
because of the diagonal form of the metric with respect to the chosen coordinates, each angle has 
value — ^. The first two terms are well defined and tend to zero when 82 tends to zero. Hence 



lim lim / kgdas = — 27r = — 27rr„. 

5,^0 52^0 igo? 

°1.°2 
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